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• We obtain the exact ground state for the Calogero-Sutherland problem in ar- 

Q\ ', bitrary dimensions. In the special case of two dimensions, we show that the 

problem is connected to the random matrix problem for complex matrices, 
provided the strength of the inverse-square interaction g = 2. In the ther- 
Qh. modynamic limit, we obtain the ground state energy and the pair-correlation 

function and show that in this case there is no long-range order. 

Calogero-Sutherland type models are a class of exactly solvable models. Starting from 
the inception till today, these models continue to be of interest, especially due to 
their exact solvability 0. Moreover, these models are related to (1 + l)-dimensional 
conformal field theory, random matrices and a host of other things || . Some variants 
of the models have attracted attention in very recent past The original model 

of Calogero described N particles in one dimension, interacting through an inverse 
square law and a two-body harmonic interaction. The latter interaction was added to 
the effect of complete discretization of the spectrum. This model gave in to an analytic 
treatment yielding an exact solution. Sutherland |J considered a variant of this 
problem where the harmonic interaction was replaced by a harmonic well containing 
the N particles. It was shown that one can exactly solve for the ground state of 
the system. The norm of the ground state wave-function was shown to be of a form 
that coincides with the joint probability density function of the eigenvalues of random 
matrices from an orthogonal, unitary or symplectic Gaussian ensemble, provided, the 
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parameter giving the strength of the inverse-square interaction was fixed to |, 1 or 2 
respectively. In an effort to generalize the original Calogero model [|l|] to dimensions 
higher than one, it was shown that some exact eigenstates including the ground 
state can be obtained for a three dimensional iV-body problem with inverse-square 
interaction, provided one also adds a long-range three-body interaction J7|, which is 
not present in one dimension. It is clearly of interest to enquire if this result of could 
be generalized to arbitrary dimensions. Perhaps a much more interesting question is 
whether one could obtain the exact ground state for the Sutherland variant of this 
model in arbitrary dimensions. One could then hope that a la one dimensional case, 
one could map the problem to some random matrix problem and obtain exact results 
for the corresponding many-body problem. The purpose of this letter is to provide 
one such example. In particular, for the special case of two dimensions, we show that 
the norm of the ground state wave-function is related to the joint probability density 
function of the eigenvalues of complex random matrices, provided, the strength of 
the inverse square interaction(g) is set to 2. Using this correspondence, we obtain 
some exact results for the ground state of a two-dimensional many-body system in 
the thermodynamic limit. In particular, in the limit N — > oo, u — > with Nuo 
fixed, we obtain the n-point correlation function. Explicit expressions for one- and 
two-point correlation functions are also given and it is pointed out that in this case 
(g = 2), unlike in one dimension, the particle density does not have the semi-circular 
form and the system does not have any long-range order. 

Let us start by describing the model in D dimensions. Following 0, we shall first 
consider the iV-body problem in D dimensions, defined by the Hamiltonian describing 
N particles interacting through an inverse square potential, a harmonic potential and 
a three body term: 

N N 1 r t- 1 N 

« = -Ev? + sS 4-+ g J2 ^ + ^ 2 £4 (i) 

i=l {i,j\i<j} V {i,j,k} ki' kj ° i=l 

i<j,i^k,j^k 

where we have chosen units such that h 2 = 2m = 1. Following 0, it is easily shown 
that some eigenstates including the ground state for the system are given by, 

V<j / \ ' b,j\Kj} J \ ' {*,j[*<j} / 

with energy 

E n = s^{2n + T D + l)uo. (3) 

Here rj is the D-dimensional position vector of the z-th particle and = — r., 
denotes the relative separation of the i-th and j-th particles, while denotes its 
magnitude. The parameters g and G are the strengths of the two- and three-body 
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potentials, respectively. Here g > — 1/2 to stop the fall to the origin. Further, 
and are two parameters determined in terms of the parameters of the Hamiltonian 

by, 



Ad = \l — 
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where LF n D denotes the Laguerre polynomial. As expected, the results are in agree- 
ment with Jl| and for D = 1 and 3 respectively. It is worth noting here that the 
case D = 2 is somewhat special, since the results are valid only for g > 0. Further, 
for any D > 1, as G and g — ► 0, the n = state corresponds to the ground state of 
the oscillator problem without the centrifugal barrier and with Bose statistics. Thus 
the situation is different from the case in one dimension ||, where, as g — ► 0, the 
ground state is that of the oscillator problem without the centrifugal barrier and with 
Fermi statistics. It is not clear whether the difference in statistics in the cases of one 
and more than one dimensions has any deeper physical significance. 

Now let us consider a Sutherland variant of the above problem, where now the 
system is described by the Hamiltonian, 

i=l {i,j\i<j} V {i,i,k} ki kj i=l 

i<j,i¥=k,j^k 

This Hamiltonian is obtained by replacing the \uj 2 
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in (|I|) with the harmonic 



well potential oj 2 J2i=i ff m the same spirit as in 0. 

The ground state eigenf unctions and eigenvalues of the system(^) are given by, 

N 



C exp 



E = [DN + N(N 
where Ad is as given in 




- r, 



(8) 



(9) 

As in the previous case, as G and g — ► 0, 
the ground state becomes that for non-interacting particles in a harmonic potential, 
without the centrifugal barrier and with Bose statistics. 
Now, rewriting \l/o from (ffl) in terms of a new variable 
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we have the probability distribution for the particles as 

*2 = c«p f-A c E% 2 )ni^-yii 2Ao > 

\ 8=1 / i<j 



(10) 



(11) 
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where C is the constant of normalization. 

In the case of D = 1, there is no three-body interaction 0, i.e. G = and Ai 
is determined from the definition^. It was shown in || that for Ai = -, 1 and 
2, the corresponding expressions for \I/q match with the distribution functions for 
the eigenvalues of random matrices from a Gaussian ensemble. In particular, they 
correspond respectively to the orthogonal, unitary and symplectic matrices ||. Using 
the well-known results derived earlier ||, Sutherland wrote down the expressions 
for density and pair-correlation functions for a many-body system in one dimension. 

Our concern for the rest of this letter will be to deal with the case of D = 2. We 
shall show that in this case also the expression for \1/q renders itself to identification 
with the distribution of eigenvalues from an ensemble of complex matrices, provided, 
9 = 2. 

In two-dimensions, the ground state energy @ becomes 

E = [2N + N(N -l)A 2 )u, (12) 
while the parameter A 2 becomes 




so that the strengths of both the interactions are equal (G = g). The form of the 
distribution (|11|) remains the same, with yj denoting two-dimensional vectors. Now, 
one can write two- vectors' as complex numbers, the two-dimensional space of positions 
of the particles becoming the complex plane, C. Let us denote the complex numbers 
corresponding to the two- vectors Vj G R 2 by Zi G C, and rewrite ([IT]) in terms of Zi 



[*o(*)f = Cexp -A 2 £ I ^ | 2 m I z i ~ z o | 2Aa , (14) 

V 1=1 / i<j 

The form (0) of is known to occur in the study of complex random matrices. 
This is precisely the form of the joint probability density function of the eigenvalues 
of matrices from an ensemble of complex matrices ||(|, provided one sets A 2 equal 
to 1. The normalization constant can be determined as in [HH, and the expression 
for \1/q with A 2 = 1 and normalized to unity is 

[*o(^)] 2 = L N ft A ex P (- E I Zi A n I ^ - Zj | 2 . (15) 

\ p=l ) \ i=l / i<j 

It is worth noting here that by fixing A 2 , one is left with one single dimensional 
variable in the problem, viz. u, which has the dimension of len ^ h 2 (recall that we are 
working in the units ft 2 = 2m = 1). Therefore, the thermodynamic limit of taking 
the area and to be large with their ratio kept finite, will be achieved by taking 
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N — ► oo, lu — ► 0, with Nu = constant. We shall show below how to take this 
singular limit. 

Following |||| one can also find out the n-point correlation functions for all n. 
Let us quote the general result. The n-point correlation function is defined as 



(N 



m\ r r N 
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(16) 



i=n+H 



For the case at hand one finds, after calculating the Van-der Monde determinant, that 
the expression for the n-point correlation function is 



7T" 



exp 



E 



det 



{i,j=l,2,-n} 

As N — > oo, the correlation functions tend to well-defined limits: 



(17) 



z n ) exp 

7T n 



E 
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In particular, the one-point correlation function, defined as, 

TZ(z) = N ^ldz 2 dz 3 ■ ■ ■ dz N , 
and interpreted as the density of the iV-particle system, is given by: 



^(C) = -exp(-C 2 )E £ T 



(18) 



(19) 



(20) 



where we have put | z \— C and omitted the suffix 1: 1Z = 1Z\. The density (|20|) is 
isotropic and does not depend on the angular coordinate 6 = arg(z). The density is 
normalized as 



2n poo 



JO 



no(d(d9 = n, 



(21) 



the total number of particles. It is interesting to note the difference between the 
density ([2t]) and the semi-circular expression for density obtained in [|J for the cor- 
responding one-dimensional system, viz. 



K(y) 



2N -y 1 



0. 



y 2 <N, 
y 2 >N. 



(22) 



Each integral in the sum in TZ(() gives a value of unity and the N terms of the sum 
adds upto N. This behavior is different from the simpler form of density in 
where one could factor out 
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yN from the density. This difference in the form of the density will be manifested 
in the different way of achieving the thermodynamic limits in the two cases. Let us 
now write down the expression for the density in terms of the original position vectors 
r. Demanding the same normalization as in (pip , this leads to the following expression 
for density 

UOr) = 4- exp (-ur 2 ) V (23) 

7TA 2 V J ^ p \ 

As noted earlier, the thermodynamic limit of the expressions is given by A" — > oo 
and uj — > 0, with Nu kept constant. However, one can see from the expression (^3]) 
that it is not meaningful to take the limit uo — > 0, since that will make the density 
vanishing. In order to take this limit meaningfully, one has to express the density in 
terms of a new variable, defined by pulling out a factor of \/N from r, 



r = VNp. (24) 
The density IZ(p) with the same normalization takes the form: 

K(p) A p( ^ v) " E «. (25) 

Denoting the density at p = by 1Z , we find that IZo = — . Thus, we see that for a 
fixed IZo, letting A^ — > oo means uo — > as ■h. 

The pair-correlation function for the system has also been obtained @, §]. It is 
given by 



TZ 2 (0 = Ki [1 - exp (-7r£ 2 Jj , (26) 

where £ = y ^^-pij. Let us note at this point that, the pair-correlation function, as 
given above in (|26"D , has apparent similarity with the pair-correlation function for the 
corresponding one dimensional case [[J, with Ai = ~, which corresponds to the case 
of orthogonal Gaussian ensemble. In both the cases the function 1Z starts off from a 
value of near £ = and grows up to unity as £ — ► oo. To facilitate discussion let 
us compare the functions Y 2 = I—IZ2 setting IZq = 1 in both the cases. The behavior 
of Y 2 for our case is 

^ 2 (0 = l-^ 2 + ^f -^ + ---, (27) 
2 o 

while for the orthogonal (Ai = ~) case of 0, one had 



y 2 (£) = 1 _ 1! + _ + . . . . (28) 
6 60 135 
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That is, while for small £, Y 2 goes in a power series with only even power of £ in our 
case, for the Ai = ^ case in one dimension, it is a power series in all powers of £. 
However, for large distances, while the Y 2 in || goes as powers of 4-, ie. 

^.(0 = ^- 1 ^ + -. (») 

with extremely mild oscillations, we have a Gaussian fall-off: Y"2 = exp (— £ 2 ), 

without any oscillations. Moreover, since the Fourier transform of a Gaussian 
function is a Gaussian only, the form factor for our system is a Gaussian, unlike that 
in §. 

For example, the 



One can also write down the higher order correlations 
three-particle correlation function in the variables p is given by: 



Ksipi, P2, Ps)/^o = 1 - exp (-tt^qpL) - ex P (-^opL) ~ ex P (-^opjL) 



+ 2 exp 



--TlTZo (Pi 2 + P23 +P31 



cos [2nn A{l, 2, 3)] 



(30) 



where .4.(1,2,3) is the area of the triangle formed by the three particles 1, 2 and 3. 
The four-body correlation function is 

7^4 (pi, P2, p 3 p 4 )/^o = 1 - exp (-tt^opL) ■ • ■ 
-7iTZ (pu + Pu)} +••• 

1 / 2 2 

,Pl2 + P23 + P31, 



+ exp 



+ 



-7T 



cos [2tt^ ^1(1,2,3)] + 



-ttKo [p\ 2 + p 2 23 + p 2 u + p\ x 



cos [2irKoA(l, 2, 3, 4)] 



(31) 



2 exp 
— 2 exp 

l z 

where ^4(1,2,3,4) = | | P13XP24 I an d ••• stands for permutations of the indices. 
It is obvious from (^) and (|31| ) as well as from the general expression ([17j) for the 
n-point correlation functions that all the distribution functions are isotropic, i.e., do 
not depend on the angular coordinate. Further, they do not show either long-range 
or quasi-long-range order and have a Gaussian fall-off at large distances. 

as that occurs in the 



The form of the distribution function (|TT|) is the same [II 



case of two-dimensional Coulomb gas |]11|,|12| , although the physics in the two cases 
are quite different. 

To conclude, in this letter, we have studied a many-body problem in two dimen- 
sions. We have written down explicit expressions for the ground state energy and the 
pair-correlation function, in case the strength of the inverse square interaction is set 
to g — 2. We then show that in this case, there is no long-range or quasi-long-range 
order. It would be very interesting if one could obtain exact results for other values 
of g and see if there is long-range order for any value of g. It will also be interesting 
if one could solve a similar many-body problem in higher dimensions, especially in 
three dimensions. But to our knowledge, the corresponding random matrix problem 
has not been worked out as yet. 
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Note Added: 

After this letter was accepted for publication, we became aware of the work of Girvin 
and MacDonald where they showed that the gauge-transformed Laughlin wave- 
function [ eq. (7) of their paper] shows off-diagonal long-range order. It then im- 
mediately follows that the Calogero-Sutherland ground state wave-function in two 



dimensions as given by (^j) [which is identical to eq. (7) of [CJ] also exhibits off- 
diagonal long-range order. 

References 

F Calogero, J. Math Phys. 10(1969)2191,2197; ibid. 12(1971)419. 
M Olshanetsky and A Perelomov, Phys. i2ep.71(1981)313; ibid. 94(1983)313. 
B Simons, P Lee and B Altshuler, Phys. Rev. Lett. 72(1994)64 . 
A Khare, J. Phys. A 29(1996)L45. 
P Ghosh, To appear. 

B Sutherland, J. Math Phys. 12(1971)246. 
F Calogero and C Marchioro, J. Math Phys. 14(1973)182. 
M L Mehta, Random Matrices, Second Edition, Academic Press, NY(1991). 
J Ginibre, J. Math Phys. 6(1965)440. 
B Jancovici, Phys. Rev. Lett. 46(1981)386. 
F Dyson, J. Math Phys. 3(1962)140,157,166. 
E Wigner, Proc 4th Can Math Cong, Toronto, 1959; p 174. 
S Girvin and A McDonald, Phys. Rev. Lett. 58(1987) 1252. 



8 



